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Abstract. This note is devoted to show how to push forward the algebraic integration 
setting in order to treat differential systems driven by a noisy input with Holder regularity 
greater than 1/4. After recalling how to treat the case of ordinary stochastic differential 
equations, we mainly focus on the case of delay equations. A careful analysis is then 
performed in order to show that a fractional Brownian motion with Hurst parameter 
■ H > 1/4 fulfills the assumptions of our abstract theorems. 
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1. Introduction 

A differential equation driven by a <i-dimensional fractional Brownian motion B = 
(B 1 , . . . , B d ) is generically written as: 

P<: y t = a+ I a(y s )dB s , te [0,T], (1) 

Oh i 

J ■ where a is an initial condition in R n , a : W 1 — > W n,d is a smooth enough function, and T is 

an arbitrary positive constant. The recent developments in rough paths analysis [41 [TBI 18] 
have allowed to solve this kind of differential equation when the Hurst parameter H of the 
fractional Brownian motion is greater than 1/4, by first giving a natural meaning to the 
integral j^a{y s ) dB s above. It should also be stressed that a great amount of information 
has been obtained about these systems, ranging from support theorems [7j to the existence 
of a density for the law of y t at a fixed instant t (see (2j [3]). 
^ . In a parallel but somewhat different direction, the algebraic integration theory (in- 

troduced in [9]), is meant as an alternative and complementary method of generalized 
integration with respect to a rough path. It relies on some more elementary and explicit 
formulae, and its main advantage is that it allows to develop rather easily an intuition 
about the way to handle differential systems beyond the diffusion case given by (DO). This 
fact is illustrated by the study of delay p[6] and Volterra [5] type equations, as well as 
an attempt to handle partial differential equations driven by a rough path [TT]. In each 
of those cases, the main underlying idea consists in changing slightly the basic structures 
allowing a generalized integration theory (discrete differential operator S, sewing map A, 
controlled processes) in order to adapt them to the context under consideration. While 
the technical details might be long and tedious, let us insist on the fact that the changes in 
the structures we have alluded to are always natural and (almost) straightforward. Some 
twisted Levy areas also enter into the game in a natural manner. 

However, all the results contained in the references mentioned above concern a fractional 
Brownian motion B with Hurst parameter H > 1/3, while the usual rough path theory 
enables to handle any H > 1/4 (see [4] for the explicit application to fBm). The current 
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paper can then be seen as a step in order to fill this gap, and we shall deal mainly with two 
kind of systems: first of all, we will show how to solve equation (pQ) when 1/4 < H < 1/3, 
thanks to the algebraic integration theory. The results we will obtain are not new, and 
the algebraic integration formalism has been extended to a much broader context in pU] 
by means of a tree-based expansion (let us mention again that the case H > 1/4 is also 
covered by the usual rough path theory) . This study is thus included here as a preliminary 
step, where the changes in the structures (new definition of a controlled path, introduction 
of a Levy volume) can be exhibited in a simple enough manner. 

Then, in a second part of the paper, we show how to adapt our formalism in order to 
deal with delay equations of the form: 

f dy t = a(y t ,y t - ri ,...,yt-r q )dB t te[0,T\, f . 

\j/t = 6, *e[-r ff ,0], [) 

where y is a R n -valued continuous process, q is a positive integer, a : M n ' 9+1 — > M. n ' d is 
a smooth enough function, B is a <i-dimensional fractional Brownian motion with Hurst 
parameter H > 1/4 and T is an arbitrary positive constant. The delay in our equation 
is represented by the family < r\ < . . . < r q < oo, and the initial condition £ is taken 
as a regular enough deterministic function on [— r q ,0]. Though this kind of system is 
implicitly considered in [12] in the usual Brownian case, and in [6] for a Hurst parameter 
H > 1/2, the rough paths techniques have only been used in this context (to the best of 
our knowledge) in p2], where a delay equation driven by a fractional Brownian motion 
with Hurst parameter H > 1/3 is considered. Our paper is thus an extension of this last 
result, and we shall obtain an existence and uniqueness theorem for equation (J2j) in the 
case H > 1/4, under reasonable regularity conditions on a and £. 

From our point of view the example of delay equations, which is interesting in its own 
right because of its potential physical applications, is also worth studying in order to see 
the kind of algebraical structures which pop out when changing the type of rough differen- 
tial system we are trying to handle. In case of a delay equation driven by a rough path of 
order 3 like ours, we shall introduce the notion of doubly delayed controlled processes, and 
have to assume a priori the existence of some doubly delayed elements of area and volume 
associated to B. This rich structure induces some cumbersome computations when one 
decides to expand all the calculations explicitly like we did. However, in the end, one also 
gets the satisfaction to see that the algebraic integration setting is flexible enough to be 
adapted naturally to many situations. Let us also mention that the infinite dimensional 
setting of [H] is avoided here, and that all our considerations only involve paths taking 
values in a finite dimensional space. 

Let us also mention that, as in other examples of fractional differential systems, an 
important part of our work consists in verifying that the fractional Brownian motion 
satisfies the assumptions of our abstract theorems. The main available tools we are 
aware of for this kind of task are based on Russo-Vallois approximations [H], analytic 
approximations of the fBm (like we did in [5]) or Malliavin calculus. We have chosen 
here to work under this latter framework, since it leads to reasonably short calculations, 
and also because it allows us to build on the previous results obtained in [16], where this 
formalism was also adopted. 

Here is how our article is structured: Section [2] is devoted to recall the basic ingredients 
of the algebraic integration setting. The diffusion case is treated at Section [3l and the 



bulk of the computations concerning delay systems can be found at Section [H Finally, 
the application to fractional Brownian motion is given at Section El 



2. Increments 

To begin with, let us present the very basic algebraic structures which will allow to 
define a pathwise integral with respect to irregular functions. 

2.1. Basic notions of algebraic integration. For an arbitrary real number T > 0, a 
vector space V and an integer k > 1 we denote by C k {V) the set of functions g : [0, T] k —> 
V, g(t 1 , . . . ,t k ) = gt!...t k such that gt x - t k — whenever ti = t i+ i for some 1 < % < k — 1. 
Such a function will be called a {k — l)-increment, and we will set C*(V) = U k >iC k (V). 
On Ck{V) we introduce the operator 5 defined as follows: 

fe+i 

5 : C k (V) -> C k+1 (V), (5g) tl ... tk+1 = ^- 1 ^' i 9t,-4,..t k+1 , (3) 

i=i 

where U means that this particular argument is omitted. A fundamental property of 5, 
which is easily verified, is that 5 o 5 = 0. We will denote ZC k (V) = CkiV) PI Ker5 and 
BC k (V) =C k (V)nim5. 

Throughout the paper we will mainly deal with actions of 5 on Cj, i = 1,2. That is, 
consider g £ C\ and h £ C 2 - Then, for any s,u, t £ [0, T], we have 

(Sg)st = 9t~ 9s, and (Sh) sut = h st - h su - h ut . (4) 

Furthermore, it is easily checked that ZC k+ i(V) = BC k (V) for any k > 1. In particular, 
we have the following property: 

Lemma 2.1. Let fc > 1 and h £ ZC k +i(V). There exists a (non unique) f £ Cjfc(V) such 
that h — Sf. 

Lemma [2711 implies that all the elements h £ C 2 (V) such that Sh = can be written 
&s h = df for some (non unique) / £ Ci(V). Thus we have a heuristic interpretation of 
S\c 2 (v) as a measure of how much a given 1-increment is far from being an exact increment 
of a function, i.e., a finite difference. 

Remark 2.1. Here is a first elementary but important link between these algebraic struc- 
tures and integration theory. Let / and g be two smooth real valued functions on [0, T] . 
Define I £ C 2 by 



(J^ dg w ^jdf v , for s,t e [0,T]. 



Then, (8I) sut = [g u - g s ][ft - fu] = (Sg) su (8f) ut . Hence we see that the operator 5 
transforms iterated integrals into products of increments, and we will be able to take 
advantage of both regularities of / and g in these products of the form SgSf. 

Let us concentrate now on the case V — R d , and notice that our future discussions 
will mainly rely on /c-increments with k < 2, for which we will use some analytical 
assumptions. Namely, we measure the size of these increments by Holder-type norms 
defined in the following way. For / £ C2(V) and n £ (0, oo), let 

\Jst\ / r \ 

SU P 77 7J7J' (°) 



s,te[o,T] 



t - sf 



and set C£(V) = {/ G C 2 (V); < 

The usual Holder spaces Cf (V) will be determined in the following way. For a continuous 
function g G Ci(V), we simply set 

IMUHIMU ( 6 ) 

where the right-hand side of this equality is defined after j5]); we will say that g G Cf (V) 
iff H^llju is finite. Notice that || • ||^ is only a semi-norm on C\(V). However we will generally 
work on spaces of the type 

C&V) = {g : [0, T] ^V;g = a, |MU < 00} , 

for a given a G V, on which \\g\\^, then becomes a norm. For h G C^{V) we set 



\\ h h, P = sup 

s,u,te[0,T] 





hsut 




M - S T £ 


— M 


P 



where the last infimum is taken over all sequences {hi G C%(V)} such that h = 
Then || • || At is easily seen to be a norm on Cs(V), and we set 

Cg(V) := {heC 3 (V); ||%<oo}. 

Eventually, let Cj + (V) = U^iC^fV), j = 1, 2, 3, and remark that the same kind of norms 
can be considered on the spaces ZC^iV), leading to the definition of some spaces ZC%(V) 
and ZCl + (V). 

With these notations in mind, the crucial point in our approach to pathwise integration 
of irregular processes is that, under mild smoothness conditions, the operator 5 can be 
inverted. This inverse is called A, and is defined in the following proposition, whose proof 
can be found in [9l [TT] : 

Proposition 2.2. There exists a unique linear map A : ZC\ + (V) — > C\ + (V) such that 

5 A = Id Z ci+(v) an d ^ = ^c\ + {v)- 

In other words, for any h G C\ + (V) such that 5h = there exists a unique g = A(h) G 
C\ + (V) such that 5g = h. Furthermore, for any \i > 1, the map A is continuous from 
ZC^iy) to C^iV) and we have 

||A%<^^INU heZC^V). 

It is worth mentioning at this point that A gives raise to a kind of generalized Young 
integral, which is a second link between the algebraic structures introduced so far and a 
theory of generalized integration: 

Corollary 2.3. For any 1-increment g G C2(V) such that 5g G C\ + , 

n 

(Id - A5)g = lim ^ g u u+1 , 

where the limit is over any partition U st = {t = s, . . . ,t n = t} of [s, t], whose mesh tends 
to zero. Thus by setting Sf = (Id — A5)g, the 1-increment Sf is the indefinite integral of 
the 1-increment g. 



We can now explain heuristically how our generalized integral will be defined. 

Remark 2.2. Let / and g be two real valued smooth functions, and define I G C 2 like in 
Remark 12 . 1L Thanks to this remark and Proposition 12.21 the following decomposition- 
recomposition for / = J df J dg holds true: 

Jdgjdf -^(6g)(6f)^fdgfdf, 

where for the second step of this construction, we have only used the fact that the product 
of increments (5g) (Sf), considered as an element of ZC 3 , is smooth enough. This simple 
procedure allows then to extend the notion of iterated integral to a non-smooth situation, 
by just applying the operator A to (5g) (5f) whenever we are allowed to do it. 

2.2. Some further notations. We summarize in this section some of the notation which 
will be used throughout the paper. 

A multilinear operator A of order /, from R dl x . . . x M. dl to M n , is denoted as an element 
A G R n > dl >-> <i i. In order to avoid tricky matrix notations, we have decided to expand all our 
computations in coordinates, and use Einstein's convention on summations over repeated 
indices. Notice that we will also use the notation A G M. dl ' d2 ' d3 ' di for a linear operator 
from M. d3,d4 to ~R dl ' d2 . We hope that this convention won't lead to any ambiguity. The 
transposed of a matrix M G M. dl < d2 is written as M* . 

For a function cp : (M n ) 9+1 — > R, we denote by <3/(^(u> , Wi, . . . ,w g ) the derivative of <p 
with respect to the i th component of Wj, for i < n and j — 0, . . . , q. 

We shall meet two kind of products of increments: first, for g G C n (M. l,d ) and h G C m (M. d ) 
we set gh for the element of C n+m ^i(MJ) defined by 

(gh)t 1 ,...,t m+n - 1 — gti,...,t n htn,~,t m+n -ll *1> • ■ • jtm+n-1 g [o,T]. (7) 

If now g G C n (M!' d ) and h G C n (R d ) we set g ■ h for the element of C n (M. 1 ) defined by 

(g ■ h) tl ,...,t n = 9t 1 ,...,t n h tl ,..,t n , h,...,t n G [0,T\. (8) 

In order to avoid ambiguities, we shall denote by M[f; Cf\ the ^-Holder norm on the 
space C j} for j = 1, 2, 3. For C G &(V), we also set AT[(; CfW] = sup < s < T |C|y 

The integral of a real valued function / with respect to another real valued function g, 
when properly defined, is written indistinctly as f fdg or J(fdg). 

3. The diffusion case 

In this section, we will recall the basic steps which allow to define rigorously and solve 
an equation of the form: 

yt = a+ [ a(y s )dx s , te [0,T], (9) 
Jo 

where a is an initial condition in W 1 , a : W 1 -> R n ' d is a smooth enough function, T is 
an arbitrary positive constant, and x is a generic (i-dimensional noisy input with Holder 
regularity 7 > 1/4. In the algebraic integration setting P, [10], this task amounts to 
perform the following steps: 

(1) Definition of an incremental operator 5 and its inverse A. 

(2) Definition of a suitable notion of controlled processes, and integration of those 
processes with respect to x. 



(3) Resolution of the equation thanks to a fixed point procedure in the space of con- 
trolled processes. 

Having dealt with the first of those points at Section I2TT1 we turn now to the second one, 
that is a definition of a useful notion of controlled processes. 

3.1. Weakly controlled processes. Before giving the formal definition of a weakly 
controlled process in the context of equation ([9]), let us recall that when the regularity 
of the noise is 7 > 1/4, the rough path setting relies on the a priori existence of an area 
(resp. volume) element x 2 (resp. x 3 ) satisfying the so-called Chen's relations: 

Hypothesis 3.1. The path R d -valued x is ^-Holder continuous with 7 > 1/4, and admits 
a Levy area and a volume element, that is two increments x 2 G '(R d,d ) and x 3 G 
C^ 1 (R d,d,d ) (which represent respectively J(dxdx) and J{dxdxdx), with the conventions 
of Section \2. 2\) satisfying: 

<5x 2 = 5x <g> 5x, i.e. (5(x 2 ) l:, ) sut = {5x l ) su {5x :, ) u t 

<5x 3 = x 2 ® 5x + 5x ® x 2 , i.e. (5(y?f k ) sut = (x*J« \5x k ) ut + (&z%(x 2 t ) J ' fc , 

for any s,u,t G [0, T], and any i,j, k G {1, . . . , d}. 

The geometrical assumption for rough paths (which is satisfied by the fractional Brow- 
nian motion in the Stratonovich setting) also states that products of increments should 
be expressed in terms of iterated integrals: 

Hypothesis 3.2. Let x 2 be the area process defined at Hypothesis \3.1l and denote by 
x 2 ' s the symmetric part of x 2 , i.e. x 2,s = |(x 2 + (x 2 )*). Then we suppose that for 
0<s<t<Twe have: 

x ?t s = \{&x)st ® (Sx) st . 

With these hypotheses in mind, the natural class of processes which will be integrated 
against x are processes whose increments can be expressed simply enough in terms of the 
increments of x: 

Definition 3.1. Let z be a process in C^iR 1 ) with k < 7 and 3k + 7 > 1, such that 
zq = a G R l . We say that z is a weakly controlled path based on x if Sz G C$(M!) can be 
decomposed into 

= ({yi6xi + (Cy ik (x 2 ) kj + r i , i.e. (Sz% = (CW^t + (C s 2 P(x 2 t ) fej r + rJ t , (10) 

for any 1 < i < I, 1 < j, k < d. Ln the previous decomposition, we further assume that 
C 1 G CI (R l ' d ) is a path with a given initial condition Co = b G R l > d , such that 5C 1 G C% (R l > d ) 
can be decomposed itself into: 

5{c iyj = {C JHx k + //'. i.e. (5(( l r) st = (O ijk (Sx k ) st + p%, 

for all s,t G [0,T], where ( 2 is a given path in Ci(R l,d,d ). Notice also that in the previous 
equations, r and p are understood as regular remainders, such that r G C| K (R') and 
p G Ci K (R l > d ). 

The space of weakly controlled paths will be denoted by Q K , a ,b(R l ), and a process z G 
QK,a,b(R l ) can be considered in fact as a triple (z,^ 1 ,^ 2 ). The natural semi-norm on 



) is given by 
N[Z] Qn,a,b 



(R l )]=Af[z;C*(R 1 )] + Af[C;C™(R 14 )} + Af[C l ;C*(R l ' d )} 
+ A/"[C 2 ; Cr(K M,d )] + N[( 2 ; C?{R l > d > d )] 
+ M[p; Ci K (R l ' d )} + Af[r; C 2 3k (R')], 



where the notations Af[g;C^(V)} and Af[(; Cj°(V)] have been introduced at Section{KM 

Remark 3.1. With respect to the case 7 > 1/3, the link between and ( 2 in the definition 
of controlled processes is new. This cascade relation between z, C 1 and ( 2 is reminiscent of 
the Heinsenberg group structure of Lyons' theory, and is really natural for computational 
purposes. 

We can now study the stability of controlled processes by composition with a regular 
function. 

3.2. Composition of controlled processes. The results of this section can be summa- 
rized into the following: 

Proposition 3.1. Assume Hypothesis \3.§\ holds true. Let z G Q K ,a,b(R l ) with decom- 
position [10}) . consider a regular function tp G Cf(R';R) and set z = <p(z), a = (p(a), 
b = di(p(a)b % . Then i G Q K ah(M)> an< ^ this latter path admits the decomposition 



Sz=(C 




(11) 



with 




f 1 = di<p(z)r l , 



r 2 = lidniM*) ■ (C 2 Y inkl ■ (C 2 ) i2i2fe2 ][(x 2 ) fc ^ ■ (x 2 ) fc ^] + \d lll2 ^)[r 11 ■ ^} 

+ [^^(^)-(cT ii -(c 2 r fe ][^-(x 2 )^] 

f 3 = 8tp(z) - d i ip(z)5z t - -dij(f{z)[8z l ■ 5z j ]. 



As far as (C Y is concerned, for 1 < j < d, it can be decomposed into 

S^y = (( 2 y k 5x k + p> 
where the remainder p can be expressed as p* = (p 1 )- 7 + (p 2 ) J , with: 



(12) 




Af[z; Q KAb (R)} < c^, T (l +M 3 [z; Q K) «,&(R Z )]) ■ 



(13) 



Proof. This proof is a matter of long and tedious Taylor expansions, and we shall omit 
most of the details. Let us just mention that we start from the relation: 

(8z) st = cp(z t ) - <f(z s ) = d i ip(z s )(5z l ) st + ^d ili2 ip(z s )(5z n ) st (5z l2 ) st 

+ <p(zt) - <fi(z s ) - dnp(z s )(5z l ) st - ^d ili2 ip(z s )(5z ll ) st (5z l2 ) s t. 



The desired decomposition (fill) is then obtained by plugging relation (fTOl) into the last 
identity, and expanding further. It should also be noticed that some cancellations occur 
due to Hypothesis 13.21 Relation (fT2]) is obtained in the same manner, and our bound (fTBl 
is a matter of standard computations once the expressions (fTTl) and (fT2l) are known. 

□ 

3.3. Integration of controlled paths. It is of course of fundamental importance for 
our purposes to be able to integrate a controlled process with respect to the driving signal 
x. This is achieved in the following proposition: 

Proposition 3.2. For a given 7 > 1/4 and k < 7, let x be a process satisfying Hypothesis 
\3.1i Let also m G Q K: b,c(M 1 ' d ) with decomposition m = b G M 1,d and 

(8m% t = (pi]) ij (5z?) st + (/<;P(x^ + rl t , l<z<d, (14) 

where p 1 G (lR d ' d ), pi = c G W. d > d , and where 5[i l G C% (IR d ' d ) can be decomposed into 

{^ l f)st = (Ar\5x k ) st + p% (15) 

with fi 2 G C?{R dAd ), p G C 2 2K (M d ' d ) ; r G C| re (R 1,d ). Define then z by z = aeR and 

5z = m i 8x i + (//Hx 2 )^ + (/i 2 )^(x 3 ) fcji + A(r i 8x i + p ij {^ 2 ) ji + 5(fi 2 ) ijk {x. 3 ) kji ) . (16) 

Finally, set 

Jirrjdx 1 ) = 5z. (17) 

Then, 

(i) z is well-defined as an element of Q K , a ,b(R), and Jirrtdx 1 ) coincides with a Rie- 
mann integral in case of some smooth processes m and x. 

(ii) The semi-norm of z in Q K) a,b(R) can be estimated as 

AT[z; Q^ b (R)] < c x ^ T {l + |% M + T^ K (\b\ Rl , d + N[m- Q kAc (M m )])}. (18) 
Furthermore, we obtain 

\\Sz\\ K < c T>x T'<- K (\b\ R ^+N[m; Q KM (R 1 ' d )]). (19) 

(iii) It holds 

Jstim'dx') 

n 

= ln lim nl EK( fal )w, +1 + Kr^UJ 1 + itij jk K,t q+ r) kji ] (20) 

llsi— >U 

q=Q 

for any < s < t < T, where the limit is taken over all partitions H st = {to — 
s, . . . , t n = t} of [s, t], as the mesh of the partition goes to zero. 



Proof. Here again, the proof is long and cumbersome, and we prefer to avoid most of the 
technical details for sake of conciseness. Let us just try to justify the second part of the 
first assertion (about Riemann integrals). 

Let us suppose then that x is a smooth function and that m G C^°(M}' d ) admits the 
decomposition flUD with p} G Cf°(R d ' d ), p 2 G C™(M. d ' d ' d ), p G C%°(R d ' d ) and r G C|°(R 1,d ). 
Then J{m l dx l ) is well-defined, and we have 

J Kdxi = m\[x\ - x\] + J [rtf u - m^dx^ 

for s < t, which can also be read as: 

Jirtfdx*) = m i 5x i + J^mW). (21) 
Let us now plug the decomposition ( TT4l ) into the expression ( |2"T1) . This yields 

J{m l dx l ) = rrjdx 1 + J {[{p}) ij Sx^dx^ + J([(n 2 ) ijk (^ 2 ) kj ]dx i ) + J^dx 1 ) 

= m'bx 1 + (^'(x 2 )* + (/i 2 p(x 3 ) fcii + Jir'dx'), (22) 

and observe that the terms m l 5x\ (/x 1 ) y '(x 2 )- ?l and (p 2 y jk (x 3 ) kjl in (1221 ) are well-defined 
provided that x, x 2 and x 3 are defined themselves. To push forward our analysis to the 
rough case, we still need to handle the term J{r l dx l ). Owing to ([22l we can write 

Jir'da?) = J[m l dx l ) - m i Sx i - {p}) ij {x 2 ) ji - (p 2 ) ijk (x s ) kj \ (23) 

and let us analyze this relation by applying 5 to both sides of the last identity. Invoking 
standard rules on the operator 5, and the fact that x satisfies Hypothesis 13.11 we end up 
with: 

5[J{r l dx 1 )} = 5{^) l] {y?Y l + 5(/i 2 p(x 3 ) fcjl - (fi 2 y ik 5x k (x 2 ) ji + r*oV, 

and thanks to the fact that S(p, 1 ) lj = (p 2 ) ljk 5x k + p lj , we obtain: 

^(rW)] = p ij {x 2 ) ji + 5(p 2 ) ijk (x 3 ) kjl + r l bx\ (24) 

Assuming now that p ij {y?) j \ 5(p?) ijk (x. 3 ) kji , r i 5x i G Cg with v > 1, then p ij (x 2 ) ji + 
5( / u 2 )^ fe (x 3 ) fe - Jl + r l 5x l becomes an element of Dom(A). Thus, applying A to both sides of 
(l24"j) and inserting the result into (|2TT) we get the expression (fT6l) of Proposition I3.2L This 
justifies the fact that (fT6j) is a natural expression for J{rri l dx % ). 

□ 

As in [15], the previous proposition has a straightforward multidimensional extension, 
which we state in the following corollary: 

Corollary 3.3. Let x be a process satisfying Hypothesis \3.1\ and let m G Q K ,b,c(M l,d ) with 
decomposition m = b G Mf' d and 

(Sm% t = (^)^(fe fc ) st +(/i 2 )^ fc2 (x 2 t ) fe2fcl +r:i ; = {£f k ^5x k *+p^ k \ (25) 

where (p, 1 )^ 1 G C?(R), (p?)^* e Cf(M), p ijkl G C 2k (R) and r ij G C 2 3k (R), /or i = 
1, . . . , I and j, ki, k 2 = 1, . . . , d Define z by z = a G ~R l and 

5z l = J(m ll dx j ) = m ij 5x j + (p 1 )^ (x 2 ) kj 

+ (^tffcifa^fcafcij + A(r*W + p ijfc (x 2 ) fcj ' + 5(/i 2 ) ijfclfc2 (x 3 ) fc2fclj ). (26) 
TTien £/ie conclusions of Proposition EH gtiZZ /io/c? m f/iis context. 



We also observe that our extended pathwise integral has a nice continuity property with 
respect to the driving path x, whose proof is also skipped here for sake of conciseness (see 
also [HI Proposition 4], and [El Proposition 3.12]). 

Proposition 3.4. Let x be a function satisfying Hypotheses \ 3.1\ and \3.Sl Suppose that 
there exists a sequence {x n ; n > 1} of piecewise C 1 -functions from [0,T] to R d such that 

lim M[x n - x- Cj(R d )} = 0, lim A/"[x 2 ' n - x 2 ; C?(R d4 )] = 0, 

n— +oo n— >oo 

and lim n _ >0O 7V[x 3 ' n - x 3 ; Cp(R dAd )} = 0. For n>l, define z n G Cf (R 1 ) in the following 
way: set z r 1 = b G R l and assume that 5z n can be decomposed into: 

S(z n Y = {C 1 ' n ) ij 5x j + ((^yjk^kj + ( r n^ S(C 1,n ) ij = (C 2 ' n Y jk Sx k + {p n ) ij , 

for 1 < i < I and 1 < j,k < d, where C 1,n G C^(R hd ) satisfies $ n = c G R M ; and 
C 2 ' n G Cf(R Wd ), p n G C 2k (R^) and r n G C 2 3k (R z ). £e£ also z be a weakly controlled 
process with decomposition ( EH) , swc/i i/iai Zo — Co = c > an ^ suppose that 

lim {J\f[z n - z;C"(R 1 )] +jV[C 1,n - ( 1 ;C™(R l ' d )] +N[( 1 > n - C 1 ;Cf(R ,,d )] 

n— »oo 

+ A/"[c 2,n - C 2 ; cr(K M,,i )] + A/"[c 2,n - C 2 ; £f (R lAd )} 

+ M[p n - p;Cl K (R 14 )] +M[r n -r-Cl K {R 1 )}} = 0. 
Finally, let up : R l -> R r ' d 6e a -function. Then 

lim N[J((p(z n )dx n ) - J{<p(z)dx); C%{R 1 ')} = 0. 

3.4. Rough diffusions equations. In this section, we shall apply the previous consid- 
erations to study differential equations driven by a rough signal, and recall that we first 
wish to solve simple equations of the form 

dy t = a(y t )dx t , y = a, (27) 

where t G [0, T], y is a R'-valued continuous process, a : R l — > R z,d is a smooth enough 
function, a; is a Revalued path and a G R z is a fixed initial condition. 

In our algebraic setting, we rephrase equation ( |27l ) as follows: we shall say that y is a 
solution to ( 1271 ). if t/o = a , y £ Qn,a,a{a) (R') and for any 0<s<t<Twe have 

(fy)* = JsMv)dx), (28) 

where the integral J(a{y)dx) has to be understood in the sense of Corollary 13.31 
With these notations in mind, our existence and uniqueness result is the following: 

Theorem 3.5. Let x be a process satisfying Hypotheses \3.1\ and \3.S\ and a : R l — > R l,d 
be a C£ -function. Then 

(i) Equation admits a unique solution y in Q K , a ,cr(a)(R') for any k < 7 such that 
3« + 7 > 1. 

(ii) The mapping (a,x, x 2 ,x 3 ) 1— > y is continuous from 

R l x Cj(R d ) x C^(R d ' d ) x C^{R dAd ) to C?(R l ), 

in the following sense: let z be the unique solution of fifty) in Q K) a,a(a)(R l ) and z 
the unique solution of fifty) in Q«,g,( 7 (a)(R i )> based on x,x, respectively. Then, there 



exists a positive constant c a>x>x depending only on a,x,x such that 

N[z - z- C?(R 1 )] < c Xii {\a -a\+ Af[x - x; CJ(R 1 )} 

+ A/"[x 2 -x 2 ;C2 7 (M d,d )] + A/"[x 3 -Sc i ;C^(R d ' d ' d )}}. 

Proof. As in (9J [TT], we first identify the solution on a small interval [0, r] as the fixed 
point of the map T : Q K: a,a(a)(R l ) —> Q.K,a,a{a) (R l ) defined by T(z) = z with z = a and 
5z = J(a(z)dx). The first step in this direction is to show that the ball 

B M = {z; z = a, Af[z; Q K , a , a( a){[0, r]; R 1 )] < M} (29) 

is invariant under T if r is small enough and M is large enough. However, due to Corollary 
13.31 and Proposition 13.11 invoking the fact that a is bounded together with its derivaties 
and assuming r < 1, we obtain 

A/TO; Q Kj a,a(a)(R 1 )) < c x {l + \a(a)\ Rl , d + r^(\a(a)\ Rl , d +Af[a(z); Q KAh (R l ' d )])} 

<c x A^ + r^Mnz)-Q KAh {R^ d )]} 

< ^{1 +t^ k (1 +M 3 [z; Q K , aMa) (R 1 )})} 

< C x , a {l + T^ K Af 3 [ Z ; Q K , a ,a(a)(R 1 )}}, (30) 

where a = a(a) and b = dia(a)a l '(a). Taking M > c Xj(T and r < r = (^4 -gg) Al, 

we obtain that c Xj(T (l + t^- k M 3 ) < M. Therefore, the ball B M defined at (EHJ) is left 
invariant by V. 

It is now a matter of standard considerations to settle a fixed point argument for Y on 
[0, r], and also to patch solutions on any interval of the form [kr, (k + l)r] for k > 1. The 
details of this procedure are left to the reader. 

□ 



4. The delay equation case 

This section is devoted to show how to change the diffusion setting in order to cover 
the case of delayed systems, having in mind to solve an equation of the form: 

/ dy t = a(y t ,yt- ri ,---,yt-r q )dxt te[0,T\, , . 

= 6, te[-r q ,0], [6i) 

where x is Revalued 7-Holder continuous function with 7 > 1/4, the function a is smooth 
enough, £ is a M n -valued 37-Holder continuous function, and < 77 < . . . < r q < 00. 
Notice that for notational convenience, we set r = and we shall also use the notation 

s{y) t = {y t . ri ,...,y t _ rq ), t e [0,T], (32) 

which means that equation (IBTT) can be written as: 

r dy t = a(y t _ ro ,s(y) t ) dx t t e [0, T], 
\y t = iu te[-r q ,0}. 

As in [16], the main ingredient in order to go from the diffusion to the delayed case will be 
the introduction of a new class of processes, namely the class of delayed controlled paths, 
which captures the structure of our equation. We shall thus first define this new class of 
paths, and see how to integrate them with respect to the driving process x. 



4.1. Delayed controlled paths. As in the diffusion case, our analysis will rely on some 
a priori increments based on our driving noise x. More specifically, we set 5(x(v)) st = 
(6x) a - v> t-v f° r s,t,v G [0,T], and we assume the following: 

Hypothesis 4.1. The path x is a ~R d -valued ■y-Holder continuous function with 7 > 
1/4, and admits two doubly delayed Levy areas and two doubly delayed volume elements. 
Namely, for v, v' G {r q , . . . , r }, we assume that there exist four paths 

x 2 (ti»,x 2 ( t /- ¥ )6C 2 2l ([0 1 T];R^), x 3 (i/, v), x 3 (i/ - v, v) G C 2 37 ([0, T]; R d ' d ' d ), 

satisfying the relations 5x 2 (V', v) = 5(x(v + v")) (g) 5(x(v)) and 

5x s (v", v) = x 2 (v", v) ® 5x + S(x(v + v")) <g> x 2 (v, r ), 

which can also be written as: 

(S(x 2 (v",v)) l3 ) sut = (5x l ) s ^ v ^ u _ v _ v n(5x 3 ) u - V) t-v 

(6{x 3 (v",v)f k ) sut = (x 2 su {v",v)f(6x k ) ut + (5x%_ v . v ^ u . v _ v ,(x 2 ut (v,r )y k , 

for v" = v' or v' — v, for any s,u,t G [0, T], and any k G {1, . . . , d}. The following 
notational simplification will also be used in the sequel: we may set x 2 (V') := x 2 (t>",w) 
whenever v = r . 

Remark 4.1. This hypothesis takes a more complex form than in [16], where the case 
7 > 1/3 was treated. However, in case of a regular process x, it should be noticed that 
the increments x 2 (f 1 ,w 2 ) and x 3 (v 1 ,v 2 ) can be defined as: 

rt—V2 rt 

x 2 st (vi,v 2 ) = / (8x(vi)) s -v2,w ® dx w , and x 3 (w 1 ,w 2 )=/ x. 2 sw (vi,v 2 ) <8 dx w , 

J S — V2 J s 

which means that x 2 (resp. x 3 ) takes the usual form of a double (resp. triple) iterated 
integral. 

As in Hypothesis 13.21 one should also express the fact that products of increments can 
be expressed in terms of iterated integrals. The following hypothesis is then easily shown 
to be a natural extension of what can be obtained in case of a smooth function x: 

Hypothesis 4.2. For v, v' G {r q , . . . , r } ; let x 2 (V, v) and x 2 (t> ' — v, v) be the area pro- 
cesses defined at Hypothesis \4-l\ Then we suppose that for all < s < t < T , we have 

x ?») = x2-v,t-vl v '> r o) and 

[Sx(v)] st ® [Sx(v')] st = y? st (v - v', v') + {x 2 st {v' - v, v))*. 

With these hypotheses in hand, the delay equation will be solved in the space of delayed 
controlled processes, which can be defined as follows: 

Definition 4.1. Let —00 < a < b < T , a given initial datum «6l" and z G ([a, b}; W 1 ) 
with k < 7 and 3/t + 7 > 1. IT^e say that z is a delayed controlled path based on x if z a = a, 
and if Sz G C 2 ([a, b\; W 1 ) can be decomposed into 

Sz l = (C^W + (C (2/) P(x 2 (ry)) fei + n\ (33) 

for all 1 < i < n, and where the index i' is summed over the set < i' < q. Just as in 
Definition \3.1l the process C 1 above has to admit the further decomposition: = [3 G W n ' d , 
where [3 has to be interpreted as another initial datum, and for 1 < j, k < d 

g^iyj = ((V> ir >y*6(x{ ril )) k + P VJ . (34) 



The regularity of the processes introduced above has to be the following: C, 1 is an element 
of C?([a,b];R n ' d ), C (2/) e C?([a,b];R n ' d ' d ), and the remainders K, p must satisfy K G 
C| K ([a, b];R n ) and p G C| K ([a, 6]; R n ' d ). 

The space of delayed controlled paths on [a,b] will be denoted by V K>a) p([a, b\; R n ), and 
a path z G T> Kja ^([a, b]; R n ) should be considered in fact as a (q + 3) -tuple (z, C\ C^ 2 '°\ ■ ■ ■ ; 
semi-norm on P KjQj/ g([a, 6]; M n ) «s £/ien gwen 6y 

JV[z;P„ )a ^([a,6];R n )] 

= Af[z; Cf([a, 6]; IT 1 )] + ATfC 1 ; ^([a, 6]; + A/'fC 1 ; CT([a, 6]; M"< d )] 
+ ^ A/"[C (2/) ; Cr([a, 6]; R nAd )] + ^ A/"[C (V) ; C?([a, 6]; R n ' d ' d )} 

i'=0 i'=0 

+ A/>; C 2 2 K ([a, b];R n > d )] + Af[TZ; Cf ([a, 6]; R")], 



where we recall that the notations J\f[g; Cf ([a, 6]; V)] and -A/"[g; C^°([a, 6]; V)] aawe 6een m- 
troduced at Section lKR 



Unfortunately, the structure above is not sufficient in order to solve the fractional delay 
equation for 7 < 1/3, and an additional notion of doubly delayed controlled processes has 
to be introduced. 

Definition 4.2. Let — 00 < a < b < T. a given initial datum a G R n and z G Cf ([a, £>]; M n ) 
w;^a k < 7 and 3k + 7 > 1. We say iAai z is a doubly delayed controlled path based on x, 
if z a = a and if Sz G C£ ([a, 6]; M n ) can 6e decomposed into 

+ (C (3/ ' J '" ) )^' fc (x 2 (r j » - r iU ,r iU )) ki + ft\ (35) 

for all 1 < i < n, 1 < j, k < d, and where the indices i',f and are summed over 

the set {1, . . . , q] and {0, 1, ... , q], respectively. As in Definition \4-l\ the processes Q 1 ' 1 ^ 
above have to admit the further decomposition: (i 1 ' 1 "* = f3^"\ where fi^"" 1 G R n,d has to be 
interpreted as another initial datum, and 



s^myj = (((w"))^( 2 ( rj ,))* + (p(o))« 

77ie regularity we ask for the processes introduced above is the following: C^- 1 ' 1 "^ is an 
element ofC?([a, b}; R n ' d ), we have (PW^W") e Cf([a, 6]; M n - d - d ), and i/ie remainders 
n, p^ r) satisfy 1Z G C 2 3K ([a, 6]; M n ) and p^") G C| K ([a, b];R n > d ). 

The space of doubly delayed controlled paths on the interval [a, b] will be denoted by 
T> K & p([a, &]; R n ), and a path z6D fi . a([a, &]; E n ) should be considered in fact as a {2q 2 + 
3g+3)-£«p/e (z, {C (1/,) }, {C (V ' j,) }, {C^"' 5 " })- ^ae natera/ semz-norm on 2^ - j([a, 6]; M n ) 



is given by 

= N[z\ Ci([a, b];R n )] + ^ A/"[C (M '' ) ; C™([a, b];R n ' d )] + ^hAf^ 1 '^; ([a, b]; R n,d )] 

i'=0 i'=0 

+ ^ Af[(( 2 > i '> jn >;C?{[a,b];R n ' d ' d )}+ ^ AT[C (2/ ' J " } ; Cf([a, 6]; R n ' d ' d )] 

i'j'=l i',j'=l 

+ ^ Af[C i3>i '' j ' ) ;Cr([a,b};R n ' d ' d )}+ Yl N[( (3 ' i '' j ' ) ;C?([a,b};R n > d > d )} 

i',j'=0 i',j'=0 

+ A/>; C 2 2K ([a, 6]; R n ' d )] + N\R,\ C 2 3K ([a, 6]; R n )]. 

As in Section^ we shall now see how delayed controlled paths behave under composition 
with a smooth map, and also how to integrate them with respect to the driving process 
x. 

4.2. Composition of delayed controlled processes. The composition of a delayed 
process with a smooth map like the function a appearing in equation (l3Tl) gives raise 
to a doubly delayed controlled process. This fact is detailed in the following proposi- 
tion, for which we recall a notation: for a function <p : (R™) 9+1 — »■ R, we denote by 
d^ip^Wo, wi, . . . , w q ) the derivative of ip with respect to the i th component of Wj, for i < n 
and j — 0, . . . , q. 

Proposition 4.1. Assume Hypothesis \4-S\ holds true. Consider < a < b < T, let a, a 

and j3, j3 be two initial conditions, respectively in R n and R n ' d , and let also ip be a function 
in C^(R n ' q+l ; R). Define a map T v on V K ^{[a,b};R n ) x V K ^{[a - r q ,b - n];R B ) by 
T v {z, z) = z, with z t = cp(z t , s(z)t) , for all t G [a, b}. Then, setting 

a = cp(a,s(z) a ), p<® = 3V(a, *(*)«)/?*, ^ = $V(«, *(*)«) (2-,.,)', 

for all 1 < %' < q, we have that T,p(z,z) G T> K& g([a,b];R). Moreover, the following cubic 
bound holds true: 

^;P^([a,6];R)] 

< tv, x>T (l + A/" 3 [z; V K>a ^{[a, b);R n )] + M 3 [z; V KAjS {[a - r„ 6 - n]; R n )]). (37) 

Proof. As in Proposition [37T], the proof of this result is based on some cumbersome Taylor 
expansions which won't be detailed here. Let us just mention briefly the decomposition 
we obtain for z: observe that 5z can be decomposed into 

(Sz)st = (2 1 ' i0 )W).-v,t-r i , + (W^.r^W^, 

+ (C^y^Urf - ry, ry))* + ti st , (38) 

where we recall that Einstein's convention on repeated indices is used, and where j, k 
are summed over the set {1, . . . , d}, and i',f over the set {0, 1, ... , q}. Furthermore, the 
expression for the coefficients £ is given by: 



and 

+ [di^(z, S (mc (2fi \r,)r%t w=n + [9^(z, S (5))(c (2j " ) )^ fc ] s % =OJ vo } , 

with the convention that C 1 ( r o) — C 1 an d C "( r °) — respectively. 

As far as the decomposition of (Q 1 ' 1 ^)^ is concerned, we obtain, for 1 < j < d: 

MC^YU = i { ^ i¥ o}(C (2/j,) ) Jfc (^) s -v-v,t-v-v 

+ (C r - I ' y"(6.r'\ ,. , + (p%>y. (39) 

For sake of conciseness, we don't include the (long) expressions we have obtained for the 
remainders TZ st , (f)g t Y here. They are also obtained via a Taylor type expansion, and the 
relation x.* t (v',v) = ^_ v t _ v (v' ,r ) (imposed by Hypothesis 14.21) turns out to be useful at 
some points of the computations. 

□ 

It should also be mentioned that, for a fixed z, the map T^-.z) : P re , Qj/ a([a, b]; R n ) — > 
T> K & z([a, b]; R) is locally Lipschitz continuous: 

Proposition 4.2. Let the notation of Proposition \4^\ prevail, and suppose that if is a 
function m C*(R n ' q+1 ; R). Let < a < b < T, letzS l \z^ £ V K ^{[a, &]; R n ) and let 
2eV K ,a,p([a-r k ,b-ri]]R n ). Then, 

M[T v {z^\ z) - T^ 2 \z);V K ^ ([a, &]; R)] 

< c^, T (l + C(* (1 \ ^ 2 \ 5)) 3 A% (1) - ^ (2) ; X>«,o,o([a, 6]; K n )], (40) 

C(^), * (2) , z) = M[z- V K ^([a -r k ,b- n]; R n )} 

+ M[z®\ V K>atf3 {[a, b);R n )] + N[z^- Z^Qa, &]; R n )] (41) 
and i/ie constant c x>V) t depends only on x, ip and T. 

4.3. Integration of delayed controlled paths. As we have seen in the previous section, 
the composition with a smooth enough function tp transforms a delayed controlled path 
into a doubly delayed controlled path. We shall see now that the integration with respect 
to x is acting in the other direction: 

Proposition 4.3. For a given 7 > 1/4 and k < 7, let x be a path satisfying Hypothesis 
4-l\ Moreover, letm £ T> K & g([a, b]; R 1,d ) such that the increments of m are given by ( fffSjJ . 



Define z by z a = a £ R and 

+ (rf W) ) y *(^(ry, r,)) kji + (£ W) )^x* fa, - ry,, ry,))^ + A st (C7), 
/or a < s < t < b, 1 < i',f < q, < < q, and where U is the increment defined by: 

u = n i 8x i + (p (< " ) ) y (x a (r i /0) i< + cKC (2/,i,) )^xV^)) fcJi 

+ rt( C -"- '" /'')^(x :5 (r ; .. - n„, n,,))^. (42) 

Finally, set 

J^rrfdx 1 ) = 8z. (43) 



Then, 

(i) J(m l dx %s ) coincides with the usual Riemann integral, whenever m and x are smooth 
functions. 

(ii) The path z is well-defined as an element ofT> K)0lj p([a, b]; R), with an initial condition 
z a = a, m a = [3 = a and a decomposition of the form 

(5z) st = + {C^F&MY* + K st , 

where the remainders 1Z G Cf K ([a, &]; E) and p G C| K ([a, 6]; M}> d ) are given by 

Kst = (Ci 2Af) y jk ^Ur^ r v )) kji + (e y ' d y j ^l(rr> - n„, r*,))** + MZ7), 
Pit = ((^y^Urf, r v )) k i + (C i -'' / ")" A '(x 2 ; (/; r - r^r^ + Kt, 
respectively. 

(Hi) The semi-norm of z in T> KjClt p(\a, b]; K) can 6e estimated as 
Af[z;V Kt0hj3 ([a,b];R)] 

< c XjTjKj7 {1 + |a| HM + (6 - a) 7 - K (|d| Ml , d + A/>; P KAj9 ([a, 6]; M M )])}. (44) 
Furthermore, the following bound also holds true: 

\\Sz\\ K < c x ,T, K , 7 (b - ay- K (\a\^ d + M[m- V KA$ ([a, b];R 1 ' d )]). (45) 

(iv) The Riemann type sums associated to our generalized integral are of the following 
form: 



MmW) = lim jh [mi kl (5x% k „ tkf+1 + (O^K^M) 



+ ($* 3 } )^l,t k ,Jr f ,r if )) kji + (Cir J T ife «„v +1 (^" - > ( 46 ) 

/or am/ a < s < t < b, where the limit is taken over any partitions H st = {to = . . . , t n = 
t} of [s,t], as the mesh of the partition goes to zero. 

Proof. As in the proof of Proposition 13.21 we shall mainly focus on the first of these 
assertions, which justifies our definition of generalized integral. Let us thus suppose for the 
moment that a; is a smooth function and that m G Ci°(M. 1 > d ) admits the decomposition (l3oT) 
with C (1/,) G C^(R d ' d ), £(2/,/) ; £(3,i'V) e C oo^ R d,d,d^ f or \ < $j < g and < < q. 
We also assume that p^ G C%°(R d ' d ) and K G C|°(R 1,d ). Then J^dx 1 ) is well-defined 
as a Riemann integral, and as for relation (1211) . one can write 

Jirrfdx 1 ) = m i 5x i + J^rrfdx 1 ). (47) 

Let us now plug the decomposition ( 1351 ) into the expression ( 1471) in order to obtain 



Jim'dx^ = m l 6x i + (C (1/ ' ) ) ii (x 2 (r i »)) J ' i + (C (W) )* ife (x 3 (7y, rv)) fci< 

+ ( C (3/'J"))yfe( x 3 (r .„ _ r .„ ; r .^)fci< + jyvdx'), (48) 

which can also be written as: 

JiTCdx 1 ) = J(Trjdx l ) - m l 5x l - (C {1/,) ) ij (x 2 (r^)) jl 

- (C (2 ' iV) P(x 3 (ry, - (C (3/V,) )^(x 3 (rj" - r^)) fc ^. (49) 



We now apply the operator 6 to both sides of the previous equation in order to get a 
suitable expression for a generalization to the rough case: 

SUCJVdx')] = Sm'Sx* + 5(C (1/ ' ) ) iJ (x 2 (r^)) ii - ((^"^[Sixin^ySx*] 

- {^^[{^{r^n^Sx' + 8{x{r v + r i 0) fc (x a (»v)) i< ] 
_ (c (3/',i"))vfc[( x 2 (r .„ _ ri „ )r . ;; ))^'^ + 5(x(r / ,)) fc (x 2 (r i »)) ii ]- 

Taking now into account the decomposition (1351 ) for m and (1361) for C^ 1 ' 1 "), we obtain: 

+ 5(C (3 ' i "' J '"))^' fc (x 3 (r i » - r>, ry/))** = 17, (50) 

where U is defined by equation ( 1421 . If we now only assume that [/ is an element of C3 
with )U > 1, one gets that [/ 6 Dom(A), which yields the decomposition (1431) for J{m l dx' L ). 
The remainder of the proof is then just made of tedious elementary estimations for the 
regularity of all the terms involved in the decomposition ([4*3l) . Finally, observe that the 
Riemann sum limit is obtained by applying Corollary 12.31 

□ 

As in Corollary 13.31 the previous Proposition can be extended easily to the multidi- 
mensional case: 

Corollary 4.4. Let m £ T> K & a([a, &]; 6e a doubly delayed controlled path, and define 
another path z by Sz z = J{m l ^dx^) for i — 1, . . . , n. Then z is well-defined as an element 
o/^K,a,/?([o, &]; R n ), and following bound holds true: 

M[z-V K ^{[a,b]-m, n )} 

< c,,t, Ki7 {1 + |a| R n.d + (6- ay- K (\a\ Rn , d + Af[m;V KA ^([a,b};R n4 )])}. (51) 

Another useful feature of our generalized integral is a continuity property with respect 
to the integrand m in J~(zdx), whose proof is omitted again for sake of conciseness: 

Proposition 4.5. Let mS l \mS 2 ^ £ V K & ^([a, b]\ IR n,a! ) 6e iwo doubly delayed controlled 
paths, and define z (1) ,z (2) £ X>«,a )( g([a, 6]; R n ) 6y 5z^ = J(m®dx), for I = 1,2. Assume 
moreover that the paths and £( 1,4 " ;2 ) m £/ie respective decompositions ( TffSj) of z^ 1 ' 

and z^ satisfy (i 1 ' 1 = Ca 1 '* ' 2 ^, and i/zai 6 — a < 1. Tnen we have: 

M[z^ -z^;V Kfifi ([a,b];R n )} 

< c x , T>K>J (b - a) KA( -^Af[mW - m ( 2 ) ; V Kfii0 ([a, b}; R n > d )}. (52) 

4.4. Rough delay differential equations. We shall now turn to the main goal of this 
section, namely the resolution of equation (I3T1) . which can be written now as: 

(8y) st = J st (a(y,s(y))dx), s,te[0,T], (53) 

with initial condition y t = £ t for t £ [— r, 0], where we recall the convention (1321) for s, and 
where the integral J has to be interpreted according to Proposition 14.31 

Before stating our main result in this direction, let us introduce a natural map, called 
T, associated to our delay equation. It is defined by V(z,z) : = z, where (recalling the 



notation T a introduced at Proposition 14.11) 52 is given as Sz = J (T a (z , z) dx) , on the 
following spaces: 

T : V K ^([a, b];R n ) x V K ^([a -r q ,b- rj; R n ) - 6]; R") 

for < a < b < T. In the previous definition, a,a,j3 and /3 stand for some initial 
conditions, with the additional compatibility condition (3 = a(a, z a _ ri , ■ ■ • > z a - r x , ^a-r 9 ), 
which shall be satisfied in our delay equation context. Notice also that, from now on, 
we shall use the convention that z s = z s = z s = £ s for s G [— r q ,0]. Since we have 
assumed £ G Cf 7 , this means in particular that, on [— r q ,0], the paths z, z, z are still 
controlled processes, whose degenerate decomposition is only given by a remainder term. 
This allows to complete easily the definition of T on intervals of the form [a, b] with a < r q . 

Remark 4.2. As in [IB], we could have handled the case of a controlled initial condition 
£. We did not consider this possibility here for sake of conciseness. 

Let us gather now some useful relations concerning the operator V we have just defined: 
first of all, by putting together inequalities (1BT1) and (137]) . one obtains 

M[V(z, 2); V K ^([a, b];M n )] < c^ tT , K „(l + Af 3 [z; V K ^([a -r q ,b- r^W 1 )}) 

x (l + (6-ar-W 3 [^P K ,^([a,6];M")]), (54) 

which means that the semi-norm of the mapping T is cubically bounded in terms of the 
semi-norm of z and z. 

Furthermore, let z (2) G V K ^([a, b];M n ) and z G T> K& ^([a - r q ,b - n]; R n ). Then, 
if b — a < 1, applying successively the inequalities ( 1521 ) and ( f40l) , we get that 

Af\T(zM, ~z) - T(zW, ~z)- V Kfi>0 ([a, b]; R n )] 

< c x ^ Ka {l + C{z^\z^ 2 \z))\b-aY h ^M^ (55) 

where C{z^\ z^ 2 \z) is defined at (|4~T]) . Therefore, for fixed z the mappings T(-,5) are 
locally Lipschitz continuous with respect to the semi-norm M\- \ P Ki o,o([ a ) % MP)]. 

With these preliminary results in hand, we can now prove our main theorem on delay 
equations: 

Theorem 4.6. Let x be a path satisfying Hypotheses \4-l\ and \4-S\ let a : R n,g+1 — > MP' d 
be a C£ -function and let £ G Cf 7 ([— r, 0]; R n ). Then 

(i) Equation (E2jJ admits a unique solution y in V K ^ o ^ r ^ )([0,T];R n ) , for any k < 7 such 
that 3k + 7 > 1 and any T > 0, where cr r (£o) — 0"(£o, £-n> ■ ■ ■ , £-rJ- 

(ii) The map 

[£, x, {x 2 (rj»)}o<i"< g , {x 3 ^', ri>)}i<i'j>< q , {x 3 (7y/ - r>, r i «)}o<i" ) j"<gJ 2/ 
is locally Lipschitz continuous from 

Cffl-r, 0];R n ) x C7([0,T];R d ) x (C 2 27 ([0, T]; R d - d )) 9+1 x (C^QO, T]; R d ' d ' d )) 2?2+29+1 

to ([0, T]; MP, in the following sense: letx be another driving rough path with correspond- 
ing delayed Levy area and doubly delayed volume element x 2 (t>) ; x 3 (t>",t>) ; respectively, 



where v" = v' or v' — v, with v,v' G {r q , . . . ,r }, and £ another initial condition. Then, 
for every M > 0, there exists a constant Km > such that the upper bound 

M[y - y; C? ([0, T};R n )} < K M [M[i - £; C^QO, T];W)\ + M[x - x; CJ([0, T];R n )} 
+ ^A/-[x 2 (r>) -x a (r i »);(?([0 i r];R* d )] 

i"=0 

+ £ ^[x»(r i /,r i O-^(r^ I ^);(? r ([0 I 7l;K** d )] 

+ ^ AT[x 3 ( V -r>,r>) -x 3 (ry, -r^,r^);C 2 37 ([0,r];M w )]|. 

i",j"=0 

/io/c?5 /or a// tuples 

(£, x, {x 2 (rj»)}o<i''< g , {x 3 (r J v, 7y)}i<i'j'< g , {x 3 0j" - r< //, r i »)}o<i"j"<?) ) 
(£,£, {x 2 (?>)}o<i"< g , {x s (rj-/,ri/)}i<i/j/< g , {x s (7y - ry/, ry)}o<i" ,/'<?)> 

satisfying the boundedness condition 
Af[£; CfQO, T];lR n )] + JV[£; ^([0, T];R n )} + A/>; #([0, T]; R")] + Af[x; Cj([0, T]; R n )} 

+ ^Ar[x 2 (r,0;C 2 27 ([0,T];M^)] + ^ AT[x 2 (r,0; C 2 27 ([0, T]; M^)] 

i"=0 i"=0 

+ ]T Ar[x 3 (r J ,,r J 0;C 2 37 ([0,T];M^)]+ £ Af[S*(r f , ry); C 2 37 ([0, T]; R^ d )] 

i',j'=l i',i'=l 

+ £ Ar[x 3 ( V -r^,^);C 2 37 ([0,T];M^)] 

i",j"=0 

+ J2 -^[x 3 ( V - ry,, r*//); C 37 ([0, T];R d ' d ' d )] < M. 

i",j"=0 

Proof. With the previous algebraic structures and computations in hand, the proof of this 
theorem follows the lines of [TU Theorem 4.2]. We shall give some hints for the proof of 
the existence-uniqueness result, which is based on a fixed point argument for the map T 
defined above, for sake of completeness. 

Without loss of generality suppose that T = Nri, where we recall that r\ is the smallest 
delay in (|53l) . We shall construct the solution of our delay equation by induction over the 
intervals [O.rJ, [0,2n], . . ., [0,iV>i]. 

Let us first show that equation (1531 ) has a solution on the interval [0, ri]. To this 
purpose, define 

/ 1 1 \t=z 

where Mi > c x = c K)7><t ,t(1 + A/" 3 [{; C^'([-r q , 0]; M n )]). In addition, choose n G [0,fi] and 
Ni G PJ such that A^n = ri, and define 

J^x^P'-l^^V!], ^' = 1,...,^. 



Finally, consider the following map: let I\i : V^^^hx,^ 1 ) -> T> K ^ 0>ar ^ o) (Ii }1 ;R n ) 
be given by z = r^i^), where 

(5z% = J at (T?(z,0dxt) 

for < s < t < T\. Notice then that if z^ 1 ' 1 ' is a fixed point of the map r^i, then z^ 1 ' 1 ' 
solves equation fl53l) on the interval Therefore, we shall prove that such a fixed point 
exists. 

First, owing to (1541 ) we get the estimate 

Af[T hl (z); T> K , M s o) (I hl] W 1 )} < d(l + rr K Af 3 [z; V K ^ o) (h,u . (56) 

Therefore, thanks to our previous choice of t±, we obtain that the ball 

B Ml = {ze T> K ,s oM s o) (h,i; K n ); M[z- V K ^ ar ^ o) (I lX , W 1 )] < M x } (57) 

is invariant under r^i. On the other hand, by changing T\ to a smaller value (and then Ni 
accordingly) if necessary, observe that also is a contraction on B AIll see (l55l) . Thus, 
applying the Fixed Point Theorem, it is easily shown that there exists a unique solution 
z^ 1 ' 1 ) to equation (l53l) on the interval 

If T\ — r±, we have thus obtained the existence and uniqueness of a solution in the 
interval [0, 7*1]. Otherwise, define the map 

r2 1 : T> (1,1) , ct,i)sU2 1; I^ n ) — ► (1.1) ^ (i.i)\(-^2 1; I^ n ) 

given by z = T 2 ,i(z), where a r (z^' 1] ) = a(z^ 1 1 ' 1 \^ , . . . ,£-r 9 _i) and 

(6?) at = J 8t {T?{z,Z)dx i ) 
for T\ < s < t < 2t\. Since T\ < r\, the following upper bound still holds true: 

and we obtain, resorting to the same fixed point argument as above, the existence of a 
unique solution z^ 2 ' 1 ^ to equation (l53l on the interval J 2 ,i- Repeating this step as often 
as necessary, which is possible since the estimates on the norms of the mappings T^, 
/ = 1, . . . , Ni are of the same type as (1541) . that is, the constant C\ in (l58|) does not change 
according to the iteration step, we obtain that z = Yli^i -^'^ 1 i s the unique solution 
to the equation (l53l) on the interval [0, ri\. 

The patching of solutions defined on different intervals of the form I^k is then a slight 
elaboration of the computations corresponding to [TBI Theorem 4.2], and this step is left 
to the reader. The continuity of the ltd map is follows also the steps of [16], except for 
the huge number of terms we have to deal with in the current situation. We prefer to 
omit this step for sake of conciseness. 

□ 

5. Application to the fractional Brownian motion 

All the previous constructions rely on the specific assumptions that we have made on 
the process x. In this section, we prove how our results can be applied to the fractional 
Brownian motion. More specifically, we first recall some basic definitions about fBm, and 
then define the delayed Levy area B 2 . We shall then turn to the definition of the volume 
B 3 , which is the main difficulty in order to go from the case H > 1/3 treated in p[6] to 
our rougher situation. 



5.1. Basic facts on fractional Brownian motion. Recall that a (i-dimensional fBm 
with Hurst parameter H G (0, 1) defined on the real line is a centered Gaussian process 



B={B t = (Bl,...,Bf); tER}, 

where B x ,...,B d are d independent 1-dimensional fBm, that is, each B % is a centered 
Gaussian process with continuous sample paths and covariance function 

R H (t, s) = E(B'iBl) = \{\t\ 2H + \s? H -\t- s\ 2H ), (59) 

for alH G {1, . . . , d}. In the sequel, all the random variables we deal with are defined on a 
complete probability space (Q, JF, P), and we assume that T is generated by the random 
variables (B t ; t G R). The fBm verifies the following two important properties: 

• Scaling property: for any c > 0, B^ = c H B./ c is a fBm, 

• Stationarity property: for any h G R, B. + h — Bh is a fBm. 

Notice that we work with a fBm indexed by R for sake of simplicity as in [16], since this 
allows some more elegant calculations for the definitions of the double delayed Levy area 
and volume, respectively. Furthermore, since the case H > 1/2 or the Brownian case 
H = 1/2 are less demanding than the rougher case, we shall mainly focus in this section 
on the range of parameter H < 1/2. 

5.1.1. Gaussian structure of B. Let us give a few facts about the Gaussian structure of 
fractional Brownian motion, following Chapter 5 of [T7]. All the considerations in this 
direction will concern a 1-dimensional fBm B, which will be enough for our applications. 

Let £ be the set of step-functions on R with values in R. Consider the Hilbert space H 
defined as the closure of £ with respect to the scalar product induced by 

%, s >]) n = R H (t', s') - R H (t', s) - R H (t, s') + R H {t, s), 

for any — oo < s < s' < +oo and — oo < t < t' < +oo, and where Rn(t, s) is given by 
( l59l) . The mapping 

%,t'] | — > B t > — B t 

can be extended to an isometry between Ti and the Gaussian space H\(B) associated with 
B. We denote this isometry by i — B(ip). 

The spaces Ti, and H^B) can be characterized more precisely in the following way: 
first, we notice that a 1-dimensional fBm defined on the real line, with H ^ 1/2, has the 
following integral representation in terms of a Wiener process W defined on R (See [2T1 
Proposition 7.2.6] for details): 

\t - s)+" 1/2 - (-s) + ~ 1/2 l dW s , t G R, (60) 



Bt -C X {H) 
where 



Cm=[f[(l + s) H ^-s H ~^] 2 ds + ±y\ (61) 

and where a + stands for the positive part of a real number a, namely a + = ]ki + (a)a. 
Using the representation (i60l ). the authors in p2| define the following stochastic integral 



of a deterministic function with respect to a 1-dimensional fBm B: 

J R (V 1 ! 2 -*/) (u)dW u , H < 1/2, 



J R (J5" 1/2 /) (u)dW u , H > 1/2, 



( X -/)W = 7^W Z-J r > (63) 



provided that the stochastic integral with respect to the Wiener process W makes sense, 
and where 

J_ r f(r) 
r(«) Ju ( r ~ u ) 

for < a < 1. The expressions ([62l and (163]) are respectively called right-sided fractional 
derivative and right-sided fractional integral on the whole real line. We remark that, in 
general, 

We also notice that 

(2?((2£/)) («) = (Z>» (22/)) («) = /(«). (64) 

When / is a function defined on an interval [a, b] with — oo<a<6<oo, extend / by 
setting /* = fl[ a ,b]- Define then 

w (-) ^ (» v) (-) = r(1 _gt^ + rd^y /; «*> 

(^/*) («) = (1«J) («) = -±-[ b j-^Jr, (66) 



T(a) A (r - «) 

for < a < 1, a < it < 6. The expressions (1651 ) and (1661) are respectively called right- 
sided fractional derivative and right-sided fractional integral on the interval [a, b] . In this 
context, as in the case of the whole line (see [20] for details and also [22]), the following 
relation holds true: 

[D_ b f) N - r(1 _ a){b _ u)a + hm r(1 _ q) / (r _ u)1+a dr. 



With these notations in hand, it is proved in [18] that the operator 

{ (v][ 2 - H f\ (u), H < 1/2, 

Cl(H) \(l H - 1/2 f)(u), H > 1/2, 
is an isometry between H and a closed subspace of L 2 (M). In fact, 

((j),i/j)n = W,/C^)l2(r), 

for all (j),tp & H. This also allows to write £>(<£>) as W(JC<p) for any <p where W(JC<p) 
has to be interpreted as a Wiener integral with respect to the Gaussian measure W. In 
particular, we have: 

E[|Z%)| 2 ] = |Mk=IIMb(R). (67) 



■5.1.2. Malliavin calculus with respect to the fBm B. Let S be the set of smooth cylindrical 
random variables of the form 



F = /(B(^), . . . , B(cp k )), <fi G H, ie{l,...,k}, 

where / G C°°{R d ' k , R) is bounded with bounded derivatives. The derivative operator 
D of a smooth cylindrical random variable of the above form is defined as the H- valued 
random variable 

i = l 1 

This operator is closable from L p (fl) into L p (fl;H). As usual, O 1,2 denotes the closure of 
the set of smooth random variables with respect to the norm 

\\F\\l 2 = E\F\ 2 + E\\DF\\ 2 n . 
In particular, considering a <i-dimensional fBm (B 1 , . . . , B d ), if D B 1 F denotes the Malli- 

12 ' 12 

avin derivative of F G with respect to B\ where ED^ denotes the corresponding 
Sobolev space, we have D B 1 B{ = #ijl(-oo,t] for i,j = l,...,d, where 5 it j denotes the 
Kronecker symbol. 

The divergence operator / is the adjoint of the derivative operator. If a random variable 
(p G L 2 (fl;H) belongs to dom(7), the domain of the divergence operator, then 7(0) is 
defined by the duality relationship 

E(FI(<p)) = E(DF,<f ) )n, (68) 

for every F G IB) 1,2 . In additon, let us recall two useful properties verified by D and /: 

• If G dom(J) and F G IB) 1 ' 2 such that Fcj) G L 2 {Vt; H), then we have the following 
integration by parts formula: 

I{F4>) = FI{4>)-(DF,(j)) H . (69) 

• If G D 1>2 (7i), D r <p G dom(J) for all r G R and {I(D r <f))} rm is an element of 
L 2 (Q;H), then 

D r I((j)) = <j) r + I(D r (j)). (70) 

One can relate the Malliavin derivatives with respect to B and W through the operator 
/C defined above. Indeed, relation (1641) shows that /C is invertible. This allows to state, 
as in the case of a 1-dimensional fBm B in an interval (see for example (TTJ Section 5.2] 
and also [1]), the following relations for the Malliavin derivative and divergence operators 
with respect to the processes B and W: 

(i) For any F G Oj^ 2 = P 1 ' 2 , we have: 

KDF = D W F, 

where D w denotes the derivative operator with respect to the process W, and Oj^ 2 
the corresponding Sobolev space. 

(ii) Dom(J) = /C _1 (Dom(J w/ )), and for any 7i-valued random variable u in Dom(J) 
we have I(u) = I w (K,u), where I w denotes the divergence operator with respect 
to the process W. 



In addition, we have D 1,2 (7Y) = (/C -1 )^ 1 ' 2 ), where L 1,2 = D 1 ' 2 (L 2 (M)), and this space 
is included in dom(I w ). Making use of the notations I w (4>) = f R 4>u dW u for any G 
dom(I w ), and I ((f)) = f R 4> u dB u for any G dom(7), we can write: 



u dB u = [ (lC<f>)(u) dW u . 
Jr 



This kind of relation also holds when one considers functions defined on an interval. 
Indeed, for some fixed — oo < a < b < oo, and H < 1/2, relation (1651 ) yields 

K dB u = [ <j> u t [aM {u) dB u = ! {K[<j>\ aM ]){u) dW u = [ (/C M 0)H dW u , 



where the operator K\ a ^ is defined by: 

(^ b] f)(u)^ T{ ^ 2 \ v]/r H f) («), for a<u<b, 

with Ci(H) defined by ( 16T1) . In case of an interval [a, 6], it should also be mentioned 
that an important subspace of integrable processes is the following: let be the set of 
step-functions on [a, b] with values in R. As in [TTJ, Subsection 5.2.3], we consider on this 
space the semi-norm 

Let Hk([o,6]) be the Hilbert space defined as the closure of £t a ' b ' with respect to the 
previous semi-norm. Then the space 7Yfc([a, b]) is continuously included in Ti, and if 
G © 1 ' 2 (H A -([a, 6])), then G Dom(J). 

5.1.3. Generalized stochastic integrals. The stochastic integrals we shall use in order to 
define our doubly delayed Levy area and volume are defined, in a natural way, by Russo- 
Vallois' symmetric approximations, that is, for a given process 0: 

J" K d°Bi = L 2 — lim 1 J' <p w (B l w+£ - B^ £ ) dw, 

provided the limit exists. It is well known that the Russo-Vallois symmetric integral 
coincides with Young's integral for H > 1/2, and with the classical Stratonovich integral 
in the Brownian case H = 1/2. Since these two cases are not very demanding from a 
technical point of view, we will focus our efforts on the case 1/4 < H < 1/2. This being 
said, for Vi G [— r, r], v 2 G [0, r], such that Vi + v 2 > 0, we will try to define the increments 
B 2 and B 3 as 

B 2 st (v u v 2 )= / d°B u ® / d°B T , i.e. (B 2 st (v u v 2 )r = / d°B{ / d°B\ 

J S — V2 J S— 1)2— VI J S—V2 J S—V2—V\ 

"t PW—V2 pu—vi 



pt rw—V2 ru—vi 

B s st ( Vl ,v 2 )= / d°B w ® / d°B u ® / d°B T , 

J S J S — V2 J S — V2—V1 

pt pW — V2 pU — Vl 

i.e. (B%( Vl ,v 2 )r k = / d°B k w / d°Bl / d°BU (71) 

J S J S — V2 J S — V2—VI 

for all i, j, k G {1, . . . ,d}, < s < t < T < 00. 

Interestingly enough, one can establish the existence of symmetric integrals thanks to 
some Malliavin calculus criterions: 



Proposition 5.1. Let <f> be a stochastic process such that 4>\ a ,b} £ ID )1 ' 2 (7i ft :([a, b])), for 
all —oo < a < b < oo. Suppose also that 

1 f b 

Tr [a:b] D(j) := L 2 - lim — / (D(p u , \ u _ £tU+e] ) H du 
s^>o ze J a 



is an 



almost surely finite random variable. Then (j) u d°B l u exists, and verifies 



! <f> u <PBl = J(# M )+lr [O)6]J D0. 

J a 



Furthermore, the following algebraic relation is trivially satisfied for this kind of inte- 
grals: 

Lemma 5.2. Let a = {a w , w E [a, b}} be a stochastic process such that its symmet- 
ric Russo-Vallois integral with respect to a 1- dimensional fractional Brownian motion B 
exists, and let F be a random variable. Then Fa is integrable with respect to B in the 
Russo-Vallois symmetric integral sense and J Fa w d°B w = F J a a w d°B w . 

We are now ready to show the existence of delayed areas and volumes with respect to 
ffim. 

5.2. Delayed Levy areas. Before we turn to statements involving increments as func- 
tions of two parameters, let us deal first with fixed times s, t: 

Proposition 5.3. Let B be a d-dimensional fractional Brownian motion, with Hurst 
parameter H > 1/4. Then, for s,t E [0, T], v\ E [— r, r], v 2 E [0, r], such that vi + v 2 > 0, 
the doubly delayed Levy area, denoted by B 2 t (i>i, i^) and defined by ( T77]) . is well defined. 
In addition, we have E[|B^(wi, v 2 )| 2 ] < c\t — s\ 4H for a strictly positive constant c = 
c h,v!,t independent of v 2 , exhibiting the following discontinuity phenomenon: we have 
lim^-j.o ch,vi,t = oo, but ch,o,t is finnite. 

Remark 5.1. The discontinuity result on ch iVi ,t alluded to above is not a surprise, and 
had already been observed in [T6] . 



Proof. As mentioned before, the case H > 1/2 is rather easy to handle, and we thus focus 
on 1/4 < H < 1/2. It should also be mentioned that Levy areas can be constructed 
in a similar way to PS], though an extra attention has to be paid in order to treat 
irregular cases, when H approaches 1/4. As a last preliminary remark, observe that, due 
to the stationarity property of the fBm we shall work without loss of generality on the 
interval [0, t — s] instead of [s — v 2 ,t — v 2 ] in the sequel, that is, B 2 t (t>i, t> 2 ) behaves as 
B^.(« 1 )=Bg it _.(« 1) 0). 

1) Casei = j andv\ > 0. Consider the process = (BJ_ V — S!_^ i )l[ ,t-s]( - )- When v\ > 0, 
the arguments in [HI Proposition 5.2] for 1/3 < H < 1/2 also hold for 1/4 < H < 1/3. 
Thus 

(B V>i)r = I B \<t>) + Tr^-^V, 
where I B * (<f>) denotes the divergence integral of <f> with respect to B l and 

T nBU-S W-s? H i if vi = 0, 

lr [0 , t _ a] iJ <p-{_ Hv 2H-l {t _ s) + l^ t _ s + Vi) 2H_ v 2H^ ifui>0 . 



In addition, one can also prove, as in [16], that 

E|(B V>i))T <c H , Vl \t-sr 



for any V\ G (0, r], where lim^^o c h,vi — °°- On the other hand, the computations above 
also show that E\{B^ t _ a {v 1} 0)) u \ 2 < c HtVl \t - s\ 4H . 

2) Case i = j and v\ < 0. When V\ < 0, we will show that 

(Bj t _>i)f = I B \<P) + Tr [0 ,-^ B >, (72) 

where now 

Tr [0 ,_ s] /J B > = H{- Vl ) 2H -\t - 8 ) + ±(\t-8 + v 1 \* H - (-^i) 2 "). (73) 

Indeed, notice that D^'(f> u = H[_ Wl)tl _ Vl ](r)l[o ) f- a ](it) and furthermore, for u € [0,t — s] and 
£ G [0, — vi], one can write 

\-^[— vi, u— vi] i e,u+e]) H 

= - (| - ui + e| 2f/ - | - «i - e| 2H + | - «! - u - £| 2H - | - - u + e| 2H ) 

= ~((-^i + e) 2 * " - ^ + I - v 1 - u - e\ 2H - | - vt - u + e\ 2H ) . 

Performing now a Taylor expansion in a neighbourhood of e = 0, we get 

(-«! + e) 2H - (-i* - e) 2H = 4H(- Vl ) 2H - l e + o {e 2 ) . 

Thus, applying the dominated convergence theorem (details are left to the reader) we 
obtain 

lim / S ±- {{-v x + e) 2H - {- Vl - e) 2H )du = H^v^-^t - s). (74) 

** 

Along the same lines, by separating the cases —V\ > t — s, < u < —V\ < t — s and 
— Vi < u < t — s, it can also be proved that 

li m I*"* ±(\- Vl - u -e\2H_\_ Vl - u + £ \ 2H )du = -{\t-s + Vl \ 2H - ( -v{) 2H ). (75) 



We now obtain (1731 ) by putting together (1741) and (175D . 

Let us bound now Tr[ t _ s ^D B1 (p from expression (l7"3l) : in the case —V\ >t — s, invoking 
the fact that, for < p < 1 and a > b > 0, the inequality a p — b p < (a — b) p holds true, 
we obtain 

|Tr [0 ,_ s] /J B >| = H{- Vl ) 2H -\t - s) + \{{-v 1 ) 2H - {- Vl -{t- s)) 2H ) 

< H{t - s) 2H + ~ {(- Vl ) 2H - {{-v x f H -{t- s) 2H )) <{t- s) 2H , 
and in the case —v\ < t — s, we also have 

1 

2 



Tr^Z^VI = H{- Vl ) 2H ~\t - s ) + -\(t-s + Vl ) 2H - (- Vl ) 21 



< H{- Vl ) 2H - l T l - 2H {t - s) 2H + \{{t- s) 2H + {- Vl ) 2H + {- Vl ) 2H ) 

< ( H (_ Vi yH-l T l-2H + ^\(t- S) 2H . 



Thus, we have found 



|Tr [0 , t „ s] /J B >| < (H^f^T 1 - 211 + fj{t- s) 2H , for all Vl G [-r,0). (76) 



We proceed now to bound the term I B \4>) in (172]) : owing to (ITU]) , we have 

Df J fi! (0) = (5^ - Bt v Jl [0 ^ s] (r) + J B * ft^.^ll^O) 

= (B l r _ Vl -BlJl [0iH (r) + / B '(ll {Bl+r , t _ a] (.))ll { - f , 1 ,t—« 1 ](r) 

= (S^ - Bljll [(¥ _ s] (r) + (5*_ s - fl^l^.^fr). (77) 

Hence, thanks to (l77j ) and using the same arguments as in the proof of [EH Propositon 
5.2], we obtain 

W B \<P)? < c H \t-s\ m , (78) 

with a constant ch > independent of v±. 

Finally, ((76]) and ([78]) imply E| (Bg it _ s (ui))"| 2 < c^Jt - s| 4H for any u x G [-r,0), and 
thus, according to our stationarity argument: 

E|(B2( Wl ,t; 2 )f| 2 < CHiW |i- s | 4 ^, (79) 

for any V\ G [— r, 0) and t> 2 G [0,r]. 

Case z ^ j. This case can be treated similarly to [TBI Proposition 5.2], and yields the 
same kind of inequality as in equation (179]) . 

Our claim E[|B^(t>i, v 2 )\ 2 } < c\t — s\ 4H now stems easily from the inequalities we have 
obtained for the 3 cases i = j and v i > 0, i = j and V\ < 0, and i 7^ j. 

□ 

We can go one step further, and state a result concerning B 2 as an increment. 
Proposition 5.4. Let B 2 be the increment defined at Proposition ^.^. Then B 2 satisfi 



e.s 



Hypothesis 4-1 and 4-2 



Proof. First, we have to ensure the almost sure existence of B 2 (v 1, v 2 ) for all s, t G [0, T}. 
This can be done by noticing that H^ t (vi,V2) is a random variable in the second chaos of 
the fractional Brownian motion B, on which all L p -norms are equivalent for p > 1. Hence 
we can write: 

E|(B 2 {vt,v 2 )y j \ p < c H , VliP \t - s\ 2pH , (80) 

for any i,j G {l,...,d} and p > 2. With the same kind of calculations, one can also 
obtain the inequality 

E|(B 2 2t > l7 v 2 )f ~ (Bl h (v u v 2 )r\ P < c Hm , p (\t 2 - t^ H + \s 2 - Sl n . 

Then, a standard application of Kolmogorov's criterion yields the almost sure definition 
of the whole family {B^ t (vi,v 2 ); s,t G [0, T]}, and its continuity as a function of s and t. 

Moreover, a direct application of Lemma 15.21 gives 

6B 2 { Vl , v 2 ) = 5(B(v 2 + Vl )) ® 5(B(v 2 )), (81) 

and Fubini's theorem for Stratonovich integrals with respect to B also yield easily Hypoth- 
esis H21 Finally, it is readily checked that B 2 (vi,v 2 ) G C^iW 14 ) for any 1/4 < 7 < H, 
Vx G [— r, r] (separating the case v\ = 0) and v 2 G [0, r]. Indeed, it is sufficient to apply 
Corollary 4 in [9] (see also inequality (90) in [16]), having in mind the bound (1801 and 
expression ( T8TT) . 

□ 



5.3. Delayed volumes. We study now the term H 3 (v 1 , u 2 ), starting from a similar state- 
ment as in Proposition I5.3t 



Proposition 5.5. Let B be a d- dimensional fractional Brownian motion, with Hurst 
parameter H > 1/4. Then, for s,t E [0, T], V\ E [— r, r], v 2 E [0, r], such that Vi+v 2 > 0, 
the doubly delayed volume, denoted by B 3 (t> 1 ,f 2 ) and defined by f7l\) . is well defined. In 
addition, we have E[|B 3 (f!, | 2 ] < c\t — s\ GH for a strictly positive constant c = Ch,t.vi,v 2 
such that it goes to oo if vi — > or v 2 — > 0, but is also well defined if v± = v 2 = 0. 

Proof. Here again, we focus on the case 1/4 < H < 1/2, and due tho the stationarity 
property of the fBm, we shall work without loss of generality on the interval [0, t — s] 
instead of [s, t] in the sequel. For notational sake, we will also set r = t — s in the 
remainder of the proof. 

1) Case i = j = k. Consider the process ip = (Bjj .(v\, v 2 ) u ) H[ 0)r ](-). We will define 
(B 3 (v 1 ,v 2 )) Ui as £ip u d Bi, which amounts to show that ip E © 1,2 (7Y/<([0, T})) and to 
compute the trace of the process ip. 

With this aim in mind, let us first compute the Malliavin derivative of ip: it is easily 
seen that 

D?^ u (82) 

= (B'r-V! ~ B -v 2 -vi)\-V2,u-v 2 ]{r)\o, T }(u) +/ B '(]l[_ t)2 _ J , 1 ,_ J ; 1 ](r)ll [ _ t)2 „ ] (-))]l[0,r](M) 
= ( B% r- Vl - B -V2-vx)\-V2,u-V2\{ r )\Q,T]{u) + (B l u _ V2 - Bl +Vi )\_ V2 _ vltU _ V2 _ Vl] (r)\ 0>r] {u). 

From this identity, one can check that ip E D 1,2 (7^x([0, T])). We will now evaluate 
J Q T ip u d°B l u by separating the Skorokhod and the trace term in the symmetric integral. 

(i) Evaluation of the trace term. We start by observing that D B ^ U can also be written 



as: 



Df Ipu = I B {\+v 1 ,u-V2\{ r )\-V2-v 1 ,u-v 2 -v 1 \{-))\o,T\{u) 

+ lB ' (%-«a-wi,-»i] ( r ) \-v 2 ,u-v 2 ] (•)) \q,t\ («) • (83) 
Apply then Fubini's Theorem in order to get 

(D B ^ U , t [u _ £:U+e] ) H du 

(/ (t[ w+vuU - V2 ],t[ u -e,u+e])ndu J dB l w (84) 

-V2— VI \J W+V2+V1 J 

/T-V2 / PT \ 
(/ (^[-vs-v^w-V!] ,\u-e,u+e])H du) dB l w , (85) 
-V 2 \Jw+V2 J 



where the last two integrals have to be interpreted in the Wiener sense, and are well- 
defined according to the criterions in [18J. 



Let us evaluate the scalar product in (I5H) : for a fixed v 2 > 0, u G [w + v 2 + v 1 ,t] 
w G [— v 2 — Vi,t — v 2 — Vi] and s G [0, U2], we can write 



1 



^2 + £ 



2H 



1 9 J-f 

V 2 — £ + 10 + 1>i — W — £ 



2iT 



+ Ui - it + £| 2H ) 



2 

= 0(^2 - ^) 2H - («2 + e) 2 * + (U - W - U! + £) 2 ^ - ( U - u; - - £) 2H ) 

= 2#( - v**" 1 + («-«;- wO 2 *" 1 ^ + o{e 2 ). 
If t>2 = 0, one can prove similarly that for e small enough, 

(l[ w+vliU _ V2 ], l[ u -e,u+e])n = 2H(U -W- Vi) 2H ~ l E + o{e 2 ) . 

This yields easily the relation 



lhn — (Ifto+fi^—ua]) %«-e,u+e]/W 



2H-1 



+ (« - w - wi) 2 ^" 1 ) if u 2 > 0, 



if (u — w — V\ 



if w 2 = 0. 



The same kind of elementary arguments work for the scalar product in expression (I85jl . 
and one obtains: 



lim ( lr_ l , 2 _ t): . m—i}-{\ 1 



l [u . £tU+£] ) n = H(-(u-w + v x ) 2H ~ l + (v 2 + Vl + u) 2H ~ l ) . 



Thus, by an application of the dominated convergence theorem (whose details are left to 
the reader) we get, for a fixed v 2 > 0, 



Tr [0 , T]J D* V 



(86) 



«2— Wl 
T— «2 



''2 



2H-1 



(r - w — v 2 — vi) + -[(t - w — Vi) 



2H „,2H- 



dBl 



(v 2 + Vl ) 2H -(r-w + Vl Y H + (r + v 2 + Wl ) 2W - (2v 2 + v l + w) 2H )dB l w , 



and for v 2 = 0, we end up with: 



T — 1)1 



(t - w - v x ) 2H dB. 



+ 



jT (v\ H -(r-w + Vl ) 2H + (r + Vl ) 2H - ( Vl + ^) 2H )rf^ 



For the remainder of the paper, the relation a < b stands for a < Cb with a universal 
constant C. Starting from equation fl86l) . let us evaluate Tr^^D B 'ip for v 2 > 0. Observe 
first that one can write E[|Tr[ , r ]I? s V| 2 ] ^ Xw=i Jh where J\ can be decomposed itself as 

j, = e[| r T F ? H^| 2 ], with ' 

= (r - w), F 2 (w) = (t + v 2 + wi) 2H - (v 2 + v x + w) 2 ^ 
F 3 (to) = (u 2 + ^i) 2H -(r-w + v 2 + Vl ) 2H , F 4 (w) = (t-w + v 2 ) 2H - v 2 2 H . 
Thus, thanks to relation (1671). we obtain: 



MIw([0,t]) - c # 



L2([0,r]) 



Furthermore, each Fi is a power function, whose fractional derivative V]/^ H Fi can be 
computed explicitly. It is then easily shown that E|Tr[ , r ]-D B V| 2 < ch,v 2 ,tt 6H , where 
ch,v 2 ,t = ch,tv^ 2H + ch- Analogously, for v 2 = 0, we get E|Tr[o, r ].D' B V| 2 < cht 6H . 

(ii) Evaluation of the Skorokhod term. We shall prove that E\I Bl (t(:)\ 2 < c HviT t 6H 
and to this aim, let us decompose ip into its Skorokhod and trace part. This gives 
W B \4>)? < 2E\I B \4> 1 )\ 2 + 2E|/^ 2 )| 2 , where 

"W—V2 



<-'2 



B l 



with 



= (BL V1 - Bt 



-V2,W—V2 



(•)• 



ip 2 (w) 
The proof that 

E|/ Bl (^)| 2 <c^ liT r 6 ^, 

where ch, Vi ,t — * oo if V\ — > but is also well defined if V\ = 0, can be obtained using the 
same arguments as for Step (i), and we then concentrate on the Skorokhod term I B \ipi). 

To estimate E|/ B °(V>i) | 2 , we use first identity fl68l) . which can be read here as E\I B * (V'i)l 2 
= E[(^i, D Bl I B * Taking into account relation (1701) . the expression (l82l) we have 
obtained for D B 'ip 1 , and the isomorphism ([671) . we end up with 

l^i)l 2 l <Qi + Q 2 + Q 3 



E 



(87) 



where Qi,Q 2 ,Q3 are respectively defined by: 

2 

Qx 



Q 2 



E 
E 

E 



^.1/2-H , 



V 



1/2-H 

~{t~V2) 



L2([0,r]) 

[ B ~m ~ B l }dB l u 



+V2 



-(t-V2-V\) 



L 2 ([-V2,T~V2]) 
2 



[B 



W — V2 



B l r , ]dBl 



+V2 + V-L 



L 2 ([— «2— «l ,r—V2—vi]) 



We now estimate those 3 terms separately, starting with Q\\ invoking the very definition 
65l) of the fractional derivative V]/ 2 H , it is easily seen that Qi < A\ + A 2 , where 



A 1 



A; 




r—V2 



[B l u - V1 ~ B\ 2 _ Vl ]dB l u 



it — r 



1-2H 



dr 



S- V2 V2 [B l u - Vl ~ B\ 2 _ Vl ]dB l u - f_ v ^\B l u _ vi - B\^ v \dB l u 

H 



(w — r) 3 / 2 



dw dr 



The term A\ is easily bounded: according to Fubini's theorem and to our previous bounds 
on B 2 , we have 



A 1 



E 



r—v 2 



[B l u _ Vl - B 



V2-V 1 \ U,JD U 



l'2 



[t — r 



,1-2H 



dr 



< c H 



AH 



[t — r 



,l~2H 



dr < ChT 



AH 



it — r 



A-2H 



dr 



°h w 
2H 



The term A 2 is a little longer to treat. However, by resorting to the same kind of tools, one 
is able to prove that A 2 < cht 6H , and gathering the estimates on A\ and A 2 , we obtain 
Qi < cht 6H as well. Finally, after some tedious computations which will be spared to the 
reader for sake of conciseness, we obtain the same kind of bound for Q 2 and Q3. 

Now one has to reverse our decomposition process: putting together our estimates on 
Qi,Q2,Q3 and plugging them into (1871) . we get E[|/' BI ('?/' 1 )| 2 ] < c H r m , with a constant 
c H > independent of vi,v 2 . Finally, gathering the bounds on the Skorokhod and the 
trace term, one obtains E[\(B 3 (v 1 ,v 2 )) iii \ 2 ] < c\t — s\ 6H . 

2) Other cases. The previous arguments and computations can be simplified to obtain 
the desired result for the case i = k ^ j and j = k 7^ i. The cases i = j 7^ k and i 7^ j 7^ k 
can be treated by means of Wiener integrals estimations. This finishes the proof of our 
claim E[\B 3 (v u v 2 )\ 2 } < c\t - s\ 6H . 

□ 

As in the case of delayed Levy areas, and with exactly the same kind of arguments, one 
can push forward the analysis in order to deal with B 3 as an increment: 

Theorem 5.6. Let B 3 be the increment defined at Proposition 15.51 Then B 3 satisfies 
Hypothesis \4-l\ Taking into account Proposition \5.j\ Theorem can thus be applied 
almost surely to the paths of the d-dimensional fBm with Hurst parameter H > 1/4. 
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